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ABSTRACT 

Splicing factor SRSF10 is known to function as a 
sequence-specific splicing activator. Here, we used 
RNA-seq coupled with bioinformatics analysis to 
identify the extensive splicing network regulated by 
SRSF10 in chicken cells. We found that SRSF10 
promoted both exon inclusion and exclusion. Motif 
analysis revealed that SRSF10 binding to cassette 
exons was associated with exon inclusion, whereas 
the binding of SRSF10 within downstream constitu- 
tive exons was associated with exon exclusion. This 
positional effect was further demonstrated by the 
mutagenesis of potential SRSF10 binding motifs in 
two minigene constructs. Functionally, many of 
SRSF10-verified alternative exons are linked to 
pathways of stress and apoptosis. Consistent with 
this observation, cells depleted of SRSF10 expres- 
sion were far more susceptible to endoplasmic re- 
ticulum stress-induced apoptosis than control cells. 
Importantly, reconstituted SRSF10 in knockout cells 
recovered wild-type splicing patterns and consider- 
ably rescued the stress-related defects. Together, 
our results provide mechanistic insight into 
SRSF10-regulated alternative splicing events in vivo 
and demonstrate that SRSF10 plays a crucial role in 
cell survival under stress conditions. 

INTRODUCTION 

Alternative splicing (AS) is a process by which different 
combinations of exons can be joined together to produce a 



diverse array of messenger RNA (mRNA) isoforms from 
a single primary transcript. This process is a common 
mechanism for controlling gene expression and generating 
protein diversity in higher eukaryotic cells. In humans, 
>95% of multiexonic protein-coding genes undergo AS 
(1); however, a precise understanding of how these 
events are regulated and coordinated in vivo is still 
lacking. Given that AS plays a key role in the regulation 
of gene expression and that aberrant splicing has been 
implicated in a wide range of human diseases (2,3), it is 
important to fully understand the regulatory mechanisms 
by which AS is regulated in vivo. 

The inclusion of an alternative exon in mature mRNA 
is largely dependent on the recognition and usage of the 
flanking splice sites by the splicing machinery. This 
involves many cis-acting RNA elements, known as 
exonic or intronic splicing enhancers and silencers, which 
are bound by different regulatory proteins. Serine/ 
arginine-rich (SR) and heterogeneous nuclear 
ribonucleoprotein (hnRNP) proteins are two well- 
characterized RNA-binding proteins that interact with 
the enhancers or silencers and subsequently stimulate or 
repress the inclusion of an alternative exon, respectively 
(4-6). However, recent genome- wide analyses revealed 
that these regulatory proteins can modulate both exon 
activation and repression in vivo, which is likely dependent 
on their binding location within the pre-mRNA (7-9). 
Therefore, it seems clear that RNA-binding proteins 
widely use this strategy to regulate splicing in the cell. 
However, splicing-sensitive arrays and CLIP-seq analysis 
of SRSF1- and SRSF2-regulated splicing revealed that 
specific splicing outcomes were not directly linked to the 
positional effect of SR protein binding in vivo. Rather, 
these outcomes were mainly determined by the balance 
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of cooperation and competition between the SR proteins 
and other regulatory proteins in vivo (10). 

SRSF10 (formally known as SRp38) is an atypical SR 
protein that functions differently from the standard SR 
proteins. Although its domain organization is similar to 
that of other SR proteins, SRSF10 was found to be unable 
to activate the splicing of a model substrate, (3-globin pre- 
mRNA, in an in vitro assay (11). SRSF10 was initially 
identified as a general splicing repressor that is activated 
by dephosphorylation (11,12). Subsequent studies showed 
that phosphorylated SRSF10 acts as a sequence-specific 
activator (13). Consistent with this role, SRSF10 was 
further demonstrated to be a key regulator of AS, 
influencing the selection of mutually exclusive exons in 
the GluR-B pre-mRNA in a sequence-specific manner 
(13,14). 

In this study, we took advantage of the previously 
generated SRSFlO-deficient (KO) chicken DT40 cells 
(12) and globally analyzed the AS events regulated by 
SRSF10 using an RNA-seq approach. We found that 
SRSF10 activates both exon inclusion and exclusion, 
and it is involved in all of the common modes of AS 
in vivo. Motif analysis revealed that SRSFlO-activated 
and -repressed exons could be distinguished by the ar- 
rangement of SRSFlO-enriched sequence motifs. 
SRSF10 binding to the alternative exon resulted in exon 
inclusion, whereas SRSF10 binding in the region of the 
downstream constitutive exon led to exon exclusion. This 
positional effect was further demonstrated by mutating 
the SRSF10 binding motifs of two minigene reporters. 
Importantly, the insertion of SRSF10 consensus sequences 
within an SRSFlO-repressed exon was sufficient to induce 
exon inclusion. Our results thus provide substantial 
insight into mechanisms of SRSFlO-regulated AS events 
in vivo. 

Functionally, a large body of literature revealed that 
many of SRSFlO-regulated alternative exons are 
associated with pathways of stress and apoptosis. 
Consistent with this finding, SRSF10 KO cells displayed 
elevated levels of apoptosis in response to endoplasmic 
reticulum (ER) stress than wild- type (WT) cells. 
Strikingly, induction of SRSF10 in KO cells recovered 
the WT splicing patterns for the majority of SRSF10- 
verified splicing events and considerably rescued the ER 
stress-related defects. These findings thus demonstrate 
that SRSF10 is required for cell survival under stress con- 
ditions, probably reflecting abilities of SRSF10 in 
regulating a large subset of stress/apoptosis-related 
splicing events. 

MATERIALS AND METHODS 

RNA sequencing 

The complementary DNA (cDNA) libraries for paired- 
end sequencing were prepared using mRNA-Seq Sample 
Prep Kit (Illumina) according to the manufacturer's 
instructions. After selecting ~200-bp fragments by 
agarose gel electrophoresis, Illumina paired-end 
sequencing adapters were ligated to the DNA fragments 
for polymerase chain reaction (PCR) amplification and 



sequencing with an Illumina HiSeq 2000 system. The 
raw sequence data and processed data have been 
submitted to Gene Expression Omnibus with accession 
number GSE53354. 

Mapping of paired-end reads 

Before read mapping, clean high-quality reads were 
obtained from the raw reads by removing the adaptor se- 
quences, reads with >5% ambiguous bases (noted as N) 
and low-quality reads containing more than half of the 
bases with qualities of <20. The clean reads were then 
aligned to the reference chicken genome (release: 
galGal4) and the corresponding recently annotated genes 
(galGal4.0) using the TopHat program (v2.0.6). In the 
alignment, preliminary experiments were performed to 
optimize the alignment parameters (-read-mismatches 3 
-read-edit-dist 4 -read-realign-edit-dist 0 -mate-inner- 
dist 20 -min-anchor 6 -coverage-search -microexon- 
search -solexal.3-quals -b2-very-sensitive) to provide the 
most information on the AS events. 

Transcriptome reconstruction 

The recent transcript annotation of the chicken genome 
sequence (galGal4.0) was obtained from ftp://ftp.ncbi. 
nlm.nih.gov/genomes/Gallus_gallus/GFF/. KO and WT 
transcriptomes were reconstructed, respectively, by 
cufflinks using the parameters (-GTF-guide 
galGal4.0.gff -frag-bias-correct galGal4.fa -multi-read- 
correct -upper-quartile-norm). Then these two transcrip- 
tomes were merged with the original version by cuffmerge 
using the parameters (-ref-gtf galGal4.0.gff -ref-sequence 
galGal4.fa KO_WT_transcriptomes.txt) to yield compre- 
hensive reannotated transcripts for subsequent AS 
analysis. 

AS analysis 

After mapping, an 'accepted_hits.bam' file was generated, 
which includes information regarding the chromosome 
position for exonic reads and exon-exon junction reads. 
All junction reads must meet two criteria: (i) at least 6 nt 
of a read must match perfectly to each of the two flanking 
regions of a potential junction site, and (ii) a junction site 
that has >3 non-redundant reads in both the WT and KO 
samples must be filtered. Afterwards, a series of Java 
programs implemented into the software named ASD 
(AS detector) were developed to fulfill the following 
tasks: (i) to reconstruct exon-clusters based on the afore- 
mentioned reannotated chicken transcriptome for identi- 
fication of five common modes of AS events for each 
exon-cluster, (ii) to count the number of junction reads 
that align either to the inclusion or exclusion isoforms in 
both the WT and KO samples, respectively, and calculate 
a P-value using junction read-counts between KO and WT 
samples by Fisher exact test, (iii) to calculate read 
coverage for the alternative exon and its corresponding 
gene in both WT and KO samples, respectively, and cal- 
culate a second P- value by Fisher exact test based on the 
alternative exon read coverage relative to its gene read 
coverage between WT and KO samples, (iv) to combine 
the two P-values to get an adjusted P-value using a 
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weighted arithmetic equation (6,15) (see Supplementary 
Method) for assessing the statistical difference of AS 
between the two samples. ASD is available at http:// 
www.novelbio.com/asd/ASD.html. 

Detection of GA-rich region and motif analysis of 
SRSFlO-regulated exons 

Based on SELEX-identified (11) and RNAcompete- 
predicted consensus sequences (16) for SRSF10, and also 
taking into account the degenerate nature of RNA- 
binding sequences, we arbitrarily used GA-rich 6-mers 
that include at least one G and two As with GA content 
>80% as potential SRSF10 binding sites. Significantly, 
enrichments of GA-rich regions were identified by 
counting the occurrences of each 6-mer (4 6 total occur- 
rences = 4096 combinations) within five segments of 
each transcript with a SRSFlO-regulated cassette exon 
(including the cassette exon and the two flanking consti- 
tutive exons and introns). The cassette exons, which were 
not affected by SRSF10, were treated as the control 
group. A hypergeometric test (Fisher's right-tailed exact 
test) was performed to detect the enrichment of each 
hexamer sequence within the sequences undergoing 
SRSFlO-regulated inclusion (or exclusion) events against 
the control sequences. Among the 4096 combinations, a 
Bonferroni-adjusted P-value (or false discovery rate, 
FDR) was calculated by multiplying the actual total 
number of tests. Subsequently, the sums of the minus 
log2-transformed P-values were calculated for GA-rich 
6-mers as a metric to determine the peaks of GA-rich 
6-mers within the five regions around the regulated 
cassette exons. Among the identified GA-rich peaks, the 
overrepresented GA-rich 6-mers (P < 0.0001) along with 
flanking 5 nt on each side (a total of 16 nt) were further 
collected and submitted to Multiple EM for Motif 
Elicitation (MEME) suite (17) for possible GA-rich 
motif analysis. The parameter setting for MEME was a 
number of single motif per sequence of 0 or 1 , a minimum 
motif width of 5, a maximum motif width of 8, a 
maximum number of different motifs of 3 and a search 
of only the given strand. Finally, MEME searching motif 
enrichments were manually edited to have a length of 6- 
mer. 

RNA isolation and RT-PCR validation 

RNA extraction, reverse-transcription and PCR amplifi- 
cation were all carried out as previously described (14). 
Briefly, PCR primers were designed to amplify two or 
multiple isoforms with different sizes. For complex 
splicing events, isoform-specific restriction digestion and/ 
or DNA sequencing was used to validate PCR product 
identity. Glyceraldehyde-3-phosphate dehydrogenase 
(GAPDH) was used as an internal control for normaliza- 
tion. Sequences for primer pairs used in this study are 
listed in the Supplementary Table S3. 

Plasmids construction, cell culture and transfection 

Minigenes were constructed by amplifying genomic 
sequences spanning exon 12-14 of KDM6A and exon 
8-10 of EIF4ENIF1, which were then inserted into 



CMV-Tag2b between the BamHl and Xhol sites, respect- 
ively. Creating SRSFlO-binding motifs on the minigene 
constructs was made as previously described (7). 
Specifically, three copies of the consensus SRSF10 
binding sequences or random sequences were inserted 
into the middle of exon 13 of KDM6A. Similar strategy 
was used to construct other plasmids used in the study. 
For deletion mutation constructs, sequences for primer 
pairs were available upon request. DT40 cells of different 
background including WT, SRSF-deficient (KO) and 
SRSFlO-deficient containing exogenously expressed hem- 
agglutinin-tagged SRSF10 (KO-re) were maintained es- 
sentially as described previously (12). Lipofectamine2000 
(Invitrogen) was used for plasmid transient transfection 
according to manufacturer's instructions. 

Cells treated with ER stress inducer, apoptosis assay and 
western blotting 

Cells were treated with the ER stress inducer (either 
thapsigargin, 0.5 uM, or tunicamycin (TM), 2ug/ml) for 
28 h. Cell apoptosis was determined by propidium iodide 
staining and flow cytometry. For recovery, after 16 h treat- 
ment with ER stress inducer, cells were collected and 
washed with phosphate-buffered saline once, then main- 
tained in fresh medium. Numbers of surviving cells were 
determined by trypan blue exclusion. Actin monoclonal 
antibodies were purchased from Santa Cruz; human 
eIF2 [pSer52] antibody was from Biosource. Western 
blotting was performed as previously described (18). 

Statistic analysis 

The data represent mean ± SD format from at least three 
independent experiments except where indicated. 
Statistical analysis was performed by Student's Mest at a 
significant level of P < 0.05. 

RESULTS 

RNA-seq and transcriptome analysis of WT and 
SRSFlO-deficient DT40 cells 

We have previously shown that SRSF10 functions as a 
sequence-specific splicing activator that is capable of 
regulating AS both in vitro and in vivo (13,14). To under- 
stand the global profile of the AS networks regulated by 
SRSF10, we used an RNA-seq approach to analyze the 
effects of SRSF10 depletion on splicing regulation in 
chicken DT40 cells (12). Briefly, WT cells and SRSF10- 
KO cells were cultured in standard medium to 90% con- 
fluence. Total RNA was isolated from these cells, and 
cDNA libraries were prepared individually by performing 
a series of procedures, including poly(A) enrichment, 
RNA fragmentation, random hexamer-primed cDNA 
synthesis, linker ligation, size selection and PCR amplifi- 
cation (Supplementary Figure SI A). The libraries were 
then sequenced in a 2 x 100-bp format using an Illumina 
HiSeq instrument. 

After quality filtering, we generated a total of 53 001 728 
and 53 138 030 clean reads for the WT and KO cells, re- 
spectively (Table 1). Then, we used TopHat (v2.0.6) to 
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Table 1. Summary of RNA-seq data mapping results 



Sample 



WT 



KO 



Clean reads 
Mapped reads 
Unique mapped reads a 
Multi-mapped reads b 
Junction reads 0 
Unmapped reads 



53 001 728 (100.0%) 53 138 030 (100.0%) 



46152 593 (87.1%) 
45 512 948 (98.6%) 

639 645 (1.4%) 
15 314481 (33.2%) 
6 847 282 (12.9%) 



45 354 730 (85.4%) 
44766437 (98.7%) 

588 293 (1.3%) 
15 505 045 (34.2%) 
7 782 321 (14.6%) 



a That include paired-end reads with an unalignable mate. 

b With multiple alignments. 

C A11 junction reads have >3 offset alignments. 



align these reads against the chicken genome sequence 
(galGal4) and its corresponding annotated genes 
(galGal4.0.gff3) (Supplementary Figure SIB). The 
paired-end reads had an approximate insert size of 
192 ± 16 bp (Supplementary Figure S2), which is consist- 
ent with the size of the cDNA selected during the library 
preparation. Given that accurate mapping of reads that 
span splice junctions is a critical component for the iden- 
tification of AS events, a series of experiments were per- 
formed to optimize the alignment parameters that would 
ensure the number of reads mapping to splice junctions 
(see the 'Materials and Methods' Section). 

After proper optimization, a total of 46 152 593 (87%) 
and 45 354 730 (85%) mapped reads were generated for 
the WT and KO cells, respectively. Almost 99% of both 
the WT and KO reads aligned to the reference genome in a 
unique manner, whereas <1% filtered as multi-mapped 
reads. Among the mapped reads, —15 314 481 WT 
(33.2%) and 15 505 045 KO (34.2%) reads were aligned 
against splice junctions (Table 1). A comprehensive anno- 
tation of the potentially novel transcripts relative to the 
known transcripts is indispensable for dissecting potential 
AS events occurring in the chicken transcriptome. To this 
end, Cufflink packages were used to reconstruct the 
chicken transcriptome based on the read distribution of 
the genomic sequences and then merge it with the current 
reference database. The final merged chicken transcrip- 
tome had a significantly expanded number of annotated 
transcripts and was sufficient for the following AS analysis 
(Supplementary Figure SIB). 

Identification of SRSFlO-affected splicing events 

There are five major modes of AS events described in 
metazoan organisms, including cassette exons, alternative 
5' splice sites, alternative 3' splice sites, mutually exclusive 
exons and retained introns (Figure 1A). To identify the 
changes in these five AS events on the loss of SRSF10, we 
developed a Java program that can search and locate sets 
of alternative exons and junctions from the genomic co- 
ordinates of all annotated chicken transcripts. Exonic 
reads and junction reads that align either to the inclusion 
or exclusion isoforms in both the WT and KO samples 
were counted, respectively. 

Comparison of the junction read counts between the 
WT and KO samples can reveal whether a change in the 
splicing pattern has occurred due to the depletion of 
SRSF10; this would be further adjusted by comparison 



A Cassette exon 
Alternative 5' splice site [ 

Alternative 3'splice site 
Mutually exclusive exon 

Retained intron [ 




Type of AS 



# of different AS events 



# of affected 



Cassette exon 


3211 


130 


Alternative 5' splice site 


1606 


16 


Alternative 3' splice site 


1738 


17 


Mutually exclusive exon 


251 


3 


Retained intron 


1089 


1 


Undefined 


2757 


6 


Total 


10652 


173 



Figure 1. Common AS events and the SRSFlO-affected splicing events 
observed in DT40 cells. (A) Diagrams of the exclusion and inclusion 
isoforms for the five modes of AS events that were examined. White 
boxes, flanking constitutive exons; black boxes, alternative spliced 
exons/regions; solid lines, splice junctions supporting the inclusion 
isoform; dotted lines, splice junctions supporting the exclusion 
isoform. (B) Summary of the different AS events identified in DT40 
cells and the splicing events affected by SRSF10 depletion, as revealed 
by analysis of RNA-seq data. 



of exon read coverage between the WT and KO cells 
(see 'Materials and Methods' Section). With an adjusted 
P- value cutoff of <0.05, we were able to identify 173 AS 
events from a total number of 10 652 splicing events that 
changed significantly in the KO samples. Specifically, the 
majority of the affected splicing events (130) belonged to 
the cassette exon category; 37 of the remaining events fell 
into other splicing modes, including 16 for alternative 5' 
splice sites, 17 for alternative 3' splice sites, 3 for mutually 
exclusive exons and 1 for intron retention. In addition, 
another six events (including coordinated two cassette 
exons and complicated splicing events) were placed into 
the 'undefined' category, as they could not be classified as 
any AS mode (Figure IB, and Supplementary Table SI). 
In summary, SRSF10 is involved in all of the common 
modes of AS, and cassette exons are the most frequent 
targets for SRSF10 regulation. 

Validation of SRSFlO-affected splicing events 

To verify whether these predicted AS events are real 
targets regulated by SRSF10, we designed primer pairs 
that were able to detect alternative exons in 45 potential 
target transcripts, which were chosen based on both 
P- values and read coverage signals. These primers were 
used to analyze the RNA isolated from the WT and KO 
cells by reverse transcriptase-polymerase chain reaction 
(RT-PCR). Although there were no or only minor 
changes in the splicing patterns in the 13 transcripts 
examined (Mest, P>0.05, date not shown), 32 were sig- 
nificantly different between the WT and KO samples 
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(Mest, P<0.05, Supplementary Table S2), indicating a 
high consistency (73%) between the RNA-seq predictions 
and RT-PCR results in our study. Unexpectedly, 24 (72%) 
events contained unannotated splice junctions. The high 
percentage of unannotated transcripts indicated that the 
currently available chicken transcript annotations are in- 
complete. Table 2 shows the top 12 splicing events that 
were dramatically affected by SRSF10 (Mest, P< 0.001). 
Figure 2 shows representative examples of these 12 events, 
with RNA-seq read coverage, RT-PCR results and quan- 
tification of their RNA products measured as the inclu- 
sion/exclusion (In/Ex) ratio. Specifically, the WT cells 
predominantly expressed cassette exon-containing 
mRNA isoforms, such as for transcripts encoding 
EIF4ENIF1 and CAST; however, in the KO cells, 
SRSF10 depletion significantly increased exon 
skipping to nearly 90% (ETF4ENIF1) and 60% (CAST) 
(Figure 2A). 

These results are in agreement with previous findings 
that SRSF10 functions as a splicing activator to promote 
exon inclusion in vivo. However, the increased inclusion of 
exon 31 of the EPRS pre-mRNA and exon 13 of the 
KDM6A pre-mRNA was observed in the SRSFlO-defi- 
cient cells (Figure 2B), suggesting that SRSF10 can also 
activate exon exclusion in vivo. In addition, the selection 
of alternative 3' and 5' splice sites was similarly affected by 
SRSF10. This was illustrated by two cases in which the 
inclusion of the BAP1 alternative 5' exon and the CDK13 
alternative 3' exon were significantly decreased on SRSF10 
depletion (Figure 2C and D). Taken together, these results 
provide strong evidence that SRSF10 can promote both 
exon inclusion and exclusion in vivo. 

Restoration of WT-splicing patterns by exogenously 
expressed SRSF10 in KO cells 

We next tested whether exogenously expressed SRSF10 
could rescue the splicing variation observed in the KO 
cells. To this end, we first examined the splicing patterns 
of 32 validated transcripts in KO-re cells, which highly 
express flu-tagged SRSF10 (12), along with those 
obtained from WT and KO cells for comparison. 



Strikingly, the induction of SRSF10 restored the WT 
splicing patterns of all of the examined transcripts with 
one exception: the inclusion of EPRS exon 11 was not 
decreased to the WT level on SRSF10 re-expression 
(data not shown). Representative RT-PCR results of the 
examined transcripts and quantification of their RNA 
products are shown in Figure 3. 

The levels of SRSFlO-activated exons were all increased 
to nearly the WT levels for the EIF4ENIF1, CAST 
and PHF20L1 transcripts on SRSF10 re-expression 
(Figure 3 A, compare lane 3 with lanes 1 and 2 for each 
target). In contrast, a decreased inclusion ratio was 
observed for the SRSFlO-repressed exons, such as for 
the transcripts encoding CEP44, L3MBTL2 and 
KDM6A (Figure 3B, compare lane 3 with lanes 1 and 2 
for each target). Similar restoration was reproduced for 
other types of alternative exons, such as the alternative 5' 
exons of SLAM and BAP1 (Figure 3C, compare lane 3 
with lanes 1 and 2 for each) and the alternative 3' exon of 
CDK13 (Figure 3D, compare lane 3 with lanes 1 and 2). 
These results strongly suggest that SRSF10 is directly 
involved in the AS regulation of these significant tran- 
scripts in vivo. 

We then determined whether the minor changes that 
were observed in several transcripts on SRSF10 deficiency 
could be rescued by the re-expression of SRSF10. 
Surprisingly, we found that these minor changes reverted 
back to the normal WT splicing levels in KO-re cells, as 
was observed for the MEF2A mutually exclusive exon 
(Figure 3E, compare lane 3 with lanes 1 and 2). This 
result provides evidence that SRSF10 also regulates the 
AS of pre-mRNAs, such as MEF2A, in vivo, although 
the output effect can be negligible. Together, these 
findings demonstrate that SRSF10 exerts differential acti- 
vation or repression effects on the various splicing events 
in vivo. 

Position-dependent activity revealed by SRSF10 binding 
motif analysis 

Next, we investigated whether the distribution of SRSF10 
binding motifs differs between SRSFlO-activated and 



Table 2. List of top 12 splicing events affected by SRSF10 



Gene symbol 


Gene description 


AS type 


Annotation of two 
spliced variants a 


EIF4ENIF1 


Eukaryotic translation initiation factor 4E nuclear import factor 1 


Exon inclusion 


A + U 


CAST 


Calpastatin 


Exon inclusion 


A + U 


CASP1 


Caspase 1, apoptosis-related cysteine peptidase 


Exon inclusion 


A + U 


PHF20L1 


PHD finger protein 20-like 1 


Exon inclusion 


A + U 


HSF2 


Heat shock transcription factor 2 


Exon inclusion 


A + A 


EPRS 


Glutamyl-prolyl-tRNA synthetase 


Exon exclusion 


A + A 


KDM6A 


Lysine (K)-specific demethylase 6A 


Exon exclusion 


A + U 


MDM4 


Mdm4 p53 binding protein homolog (mouse) 


Exon exclusion 


A + U 


CDK13 


Cyclin-dependent kinase 13 


Ait y SS 


A + U 


SLTM 


SAFB-like transcription modulator 


Alt 5' SS 


A + U 


BAP1 


BRCA1 -associated protein- 1 


Alt 5' SS 


A + U 


RBBP5 


Retinoblastoma binding protein 5 


Alt 5' SS 


A + U 



Alt; alternative: SS; splice site. 

a Based on chicken NCBI and Ensemble databases, 'U' stands for unannotated splice transcripts, whereas A stands for annotated transcripts. 
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EIF4ENIF1 




RNA-seq reads 



B 

EPRS 



RNA-seq reads 



WT 
KO 



6.0 n 
4.0 
2.0 
0.0 



WT KO 



CAST 



KDM6A 



84nt 



RNA-seq reads 
WT f| 



KO 



C 

BAP1 



RNA-seq reads 



WT 
KO 




WT KO 




WT KO 



RNA-seq reads 



WT 
KO 




RNA-seq reads 

1 



WT. 
KO 



I 



1.5 
1.2 
0.9 
■0.6 
0.3 
0 



it 



WT KO 




WT KO 



Figure 2. SRSF10 can promote both exon inclusion and exclusion in vivo. SRSFlO-activated cassette exons (A), SRSFlO-repressed cassette exons 
(B), SRSFlO-activated alternative 5' exons (C) and SRSFlO-activated alternative 3' exons (D) were identified by RNA-seq analysis and validated by 
RT-PCR. Each case is schematically diagrammed with mapped RNA-seq reads that cover the corresponding exons (left panels). Black boxes, 
alternative exons; black lines, introns; white boxes, flanking constitutive exons. The primer pairs used for RT-PCR are indicated by arrows. 
SRSFlO-regulated AS changes were analyzed by RT-PCR (right bottom panels). Quantification of the RNA products was measured as the inclu- 
sion/exclusion (In/Ex) ratio, as shown in the top right histogram. Values shown are the mean ± SD, n = 3. Significance for each detected change was 
evaluated by Student's /-test. 



-repressed cassette exons. To this end, we simply defined 
GA-rich hexamers as potential SRSF10 binding sites 
based on the SELEX-identified (11) and RNAcompete- 
predicted (16) sequences for SRSF10 (see 'Materials and 
Methods' Section). We then searched five transcript 
regions, including the cassette exon, the two flanking con- 
stitutive exons and the two flanking introns, for GA-rich 
hexamer enrichments among 78 transcripts with SRSFlO- 
activated exons and 52 transcripts with SRSFlO-repressed 
cassette exons. Another 80 SRSFlO-unaffected events 
served as controls. As shown in Figure 4 A, no GA-rich 
hexamers were found in the upstream exon or in the two 
flanking introns for both groups. However, some 
overrepresented GA-rich hexamers (e.g. GAAAGT GA 
AG A A and AG A A AG) were only found within the 
SRSFlO-activated exons, whereas other overrepresented 
GA-rich hexamers (e.g., GAAAAG AGAAA and GAA 
AAT) were found in the flanking downstream constitutive 



exons for both groups (Figure 4B). These overrepresented 
GA-rich 6-mers was further analyzed by MEME for 
possible SRSF10 binding motifs (Figure 4C; see 
'Materials and Methods' Section). 

In summary, SRSFlO-activated exons showed a pre- 
dominant enrichment of GA-rich hexamers within 
the exons themselves over the flanking downstream con- 
stitutive exon (Figure 4A). However, SRSFlO-repressed 
exons were associated with the enrichment of GA-rich 
hexamers only in the downstream constitutive exon 
(Figure 4A). Thus, the binding of SRSF10 to the alterna- 
tive exon was associated with exon inclusion, whereas the 
binding of SRSF10 within the downstream constitutive 
exon was associated with exon exclusion. This positional 
effect as to inducing exon inclusion on the cassette exon 
and skipping on the downstream exon was only previously 
demonstrated with typical SR protein SRSF1 using MS2 
tethering approach (19). 
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Figure 3. Exogenous expression of SRSF10 in KO cells restores WT splicing patterns of SRSFlO-targeted transcripts. SRSFlO-activated cassette 
exons (A), SRSFlO-repressed cassette exons (B), alternative 5' exons (C), alternative 3' exons (D) and mutually exclusive exons (E) were analyzed for 
splicing changes in WT, KO and KO-re cells by RT-PCR (left panels). The primers used for RT-PCR are the same as in Figure 3. Bar graphs show 
the ratios of inclusion versus exclusion (A-D) or exon a versus exon b (E) (right panels). Values shown are the mean ± SD, n = 3. The significance of 
each detected change was evaluated by Student's /-test. 



Mechanistic insights into SRSFlO-regulated exon 
inclusion and exclusion 

For SRSFlO-activated exons, potential SRSF10 binding 
motifs were observed predominantly within cassette exons 
over downstream exons. We thus reasoned that the 
binding of SRSF10 to cassette exons is associated with 
exon inclusion. To test this, we analyzed individual cases 
that displayed SRSFlO-dependent exon inclusion, such as 
EIF4ENIF1 exon 9. We constructed a minigene plasmid 
in which a genomic DNA fragment spanning EIF4ENIF 1 
exons 8-10 was preceded by the CMV promoter and 
followed by an SV(40) polyA site (Figure 5A, top). 
Splicing was assayed following transient expression in 



DT40 cells. Total RNA was extracted from the transfected 
cells, and exon inclusion/exclusion was analyzed by 
RT-PCR (Figure 5B). Consistent with its endogenous 
splicing pattern, the EIF4ENIF1 exon 9 was almost fully 
included in the WT cells (lane 1), whereas SRSF10 deple- 
tion significantly decreased the inclusion of exon 9 to 
nearly 30-40% (lane 2), indicating that inclusion of 
EIF4ENIF1 exon 9 was SRSF10 dependent. We identified 
four GA-rich clusters: GA1 and GA2 within exon 9 
(Figure 5A, bottom), whereas GA3 and GA4 within 
exon 10 (Supplementary Figure S3A). 

We next examined, in more detail, the role of the indi- 
vidual GA elements in exon inclusion. To this end, we 
introduced a single GA deletion into the minigene 
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Figure 4. Hexamer and MEME analysis of SRSFlO-regulated exon inclusion and exclusion. (A) The sum of the minus log2-transformed P- 
values of GA-rich hexamers within the five regions around regulated cassette exons was compared with control cassette exons that are not 
affected by SRSF10. The red line represents SRSFlO-mediated exon inclusion compared with the control (left), and the blue line represents 
SRSFlO-mediated exon exclusion compared with the control (right). There is no obvious peak in the upstream exon or the two flanking introns 
around the cassette exons. (B) A representative list of top GA-rich hexamers was shown for SRSFlO-activated cassette exons, -activated 
downstream constitutive exons and -repressed downstream constitutive exons, respectively. FC is fold change, P-value stands for statistically 
significance of GA-rich hexamers between SRSFlO-regulated group and controls and FDR indicates false discovery rate. Shown in brackets are 
the number of overrepresented GA-rich 6-mers and the percentage of exons containing at least one overrepresented 6-mer. (C) Potential GA-rich 
motifs derived from the aforementioned overrepresented GA-rich 6-mers using MEME (P < 0.0001). 



construct and analyzed the effects of each mutant on the 
exon inclusion in vivo. As shown in Figure 5C, deletion of 
GA1 led to nearly 80% exon exclusion, which is similar to 
the effects caused by SRSF10 deficiency (compare lanes 3 
and 4 with lanes 1 and 2). In contrast, the other three 
mutants displayed no effect on exon exclusion and still 
remained responsive to SRSF10 deficiency (compare 
lanes 5-6 and Supplementary Figure S3B). Taken 
together, these results provide strong evidence that the 
GA1 element has SRSFlO-dependent enhancer activity. 

To provide direct evidence that SRSFlO-dependent 
enhancer activity within the GA1 fragment can lead to 
exon inclusion, we made two other constructs in which 
three copies of the SRSF10 consensus motif or random 



sequences were inserted in the middle exon of the AGA1 
plasmid (Figure 5B). When SRSF10 binding motifs were 
present in the AGA1 plasmid, exon inclusion was restored 
nearly to the WT level, compared with the control 
sequence (Figure 5C, compare lane 9, lane 7 and lane 1). 
However, the differential effects on exon inclusion 
between the two inserted sequences were completely di- 
minished on SRSF10 depletion (Figure 5C, compare 
lane 10 and lane 8). 

For SRSFlO-repressed exons, a key insight from the 
motif analysis was the enrichment of SRSF10 binding 
motifs within the downstream constitutive exon. To inves- 
tigate the possible role of downstream SRSF10 binding 
motifs in the regulation of exon exclusion, we selected 
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Figure 5. Characteristics of SRSFlO-regulated exon inclusion and exclusion. (A) Schematic representation of the EIF4ENIF1 minigene construct 
(top). The intact genomic region of EIF4ENIF1 from exon 8-10 was cloned into the CMV-Tag2b vector. The exon 9 sequences are shown at the 
bottom, with two GA-rich motifs (GA1 and GA2) marked with a different color. The PCR primers used to detect the alternatively spliced products 
are indicated by arrows. (B) Three copies of the SRSF10 consensus binding motif or random sequences were inserted within the middle exon of the 
GA1 -deletion construct. (C) RT-PCR analysis of the EIF4ENIF1 minigene or each of the indicated mutants transfected in WT or KO cells, 
respectively. The percentage of exon inclusion is shown at the bottom of the gel. Values shown are the mean ± SD, n = 3. (D) Schematic repre- 
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or KO cells, respectively. The histogram below shows mean percentage of exon inclusion ± SD, n = 3. 



an exon from the KDM6A gene, which displayed exon 
exclusion in the WT cells but showed increased inclusion 
on SRSF10 depletion (Figure 2B). This effect was 
reproduced in the transient transfection of the WT or 
KO cells with a KDM6A minigene construct (Figure 5F, 
lanes 1 and 2). Although potential SRSF10 binding motifs 
were not observed within cassette exon 13, two GA-rich 
elements (named GA5 and GA6) were identified within 
downstream constitutive exon 14 (Figure 5D, bottom). 
Although the deletion of GA5 had little effect on exon 
inclusion, the deletion of GA6 increased exon inclusion 
from 20% to nearly 40% (Figure 5F, compare lanes 5 
and 6 with lanes 3 and 4). These results suggest that the 
GA6 motif partially contributes to exon 13 exclusion. 
Importantly, insertion of SRSF10 consensus motifs in 
the exon 14 restored the WT splicing pattern of AGA6 
substrate (Figure 5E, and F, compare lanes 7-10 with 
lanes 1-2), further demonstrating that SRSFlO-dependent 



binding motifs within the constitutive exon resulted in 
exon skipping. 

Induction of exon inclusion by the manipulation of 
SRSF10 binding motifs 

We next considered whether the introduction of SRSF10 
binding sequences in a naturally excluded exon would be 
sufficient to induce inclusion. To this end, mutations were 
made in the KDM6A minigene construct to allow for the 
insertion of three copies of the SRSF10 consensus motif or 
random sequences just downstream of the 3' splice site in 
the cassette exon (Figure 5E). As shown in Figure 5F, 
~90% of the mutated exon with the inserted SRSF10 
binding motifs was included, compared with the ~20% 
exon inclusion observed with the control sequence 
(Figure 5F, compare lane 13 and 11). Moreover, the 
effect on exon inclusion critically depended on the 
motifs that were directly bound by SRSF10, as it was 
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diminished in response to SRSF10 depletion (Figure 5F, 
lane 14). Taken together, these results suggest that the 
activation activity of SRSF10 binding sequences within 
the cassette exon is dominant over the repression activity 
of the downstream exon. 

SRSF10 protects DT40 cells from ER stress-induced 
apoptosis 

Through an extensive literature review, we learned that a 
large subset of genes have functional annotations related 
to stress and apoptosis among SRSFlO-regulated AS 
events, such as genes of HSF2 (20), CASP1 (21), CAST 
(22), SLTM (23), EIF4ENIF1 (24,25), FKBP14 (26), 
BAP1 (27), SEC16A (28) or MDM4 (29). Because 
SRSF10 is not required for cell survival (12), we ques- 
tioned whether SRSF10 favors cell growth under stress 
conditions by regulating these splicing events. To this 
end, we treated WT, KO and KO-re cells with different 
stress inducers and observed the responses of these cell 
lines to each treatment. Although there were no obvious 
differences observed among these cell lines in response to 
most stresses, ER stress induced by either thapsigargin or 
TM exerted greater deleterious effects on the KO cells 



compared with the WT and KO-re cells. As shown in 
Figure 6A, significant apoptosis was observed in the KO 
cells after 28 h of treatment with TM, as measured by 
propidium iodide staining and flow cytometry. 

We further measured the ability of these cell lines to 
recover after ER stress. Whereas both the WT and KO- 
re cells resumed rapid growth after an ~l-day lag, the KO 
cells showed considerably delayed recovery, resuming 
growth only after at least a 3-day lag (Figure 6B). These 
results reflect the loss of SRSF10 rather than an unrelated 
property of the cells because (i) another independent 
SRSF10-KO cell line behaved identically to the first 
(data not shown) and (ii) stable expression of exogenous 
SRSF10 (KO-re cells) considerably rescued the ER stress- 
related defects. In addition, in contrast to heat shock- 
induced rapid dephosphorylation (12), SRSF10 remained 
phosphorylated and dephosphorylated SRSF10 was not 
detected, regardless of short-term or long-term treatment 
of ER stress (Figure 6C). Consistent with this finding, 
none of SRSFlO-dependent exon inclusion or skipping 
was reversed under ER stress conditions (Supplementary 
Figure S4). Therefore, we concluded that it is phosphory- 
lated SRSF10 that plays a critical role in cell survival after 
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Figure 6. SRp38 protects DT40 cells from ER stress-induced apoptosis. (A) WT, KO and KO-re cells were treated with the ER stress-inducer 
tunicamycin (TM, 2ug/ml) for 28 h. Apoptosis was measured by propidium iodide staining and flow cytometry (right panel). (B) Cells were treated 
with TM (2ug/ml) for 16 h and allowed to recover in medium for 1-3 days. The number of surviving cells was determined by trypan blue exclusion. 
The error bars represent the SD. (C) WT and KO cells were subjected to TM-induced ER stress for the indicated times. Immunoblotting was 
performed with whole cell extracts using antibodies against SRSF10, phospho-eIF2A and P-actin. 



Nucleic Acids Research, 2014, Vol. 42, No. 6 4029 



ER stress, perhaps reflecting its ability to regulate many 
stress- and apoptosis-related AS events. 



DISCUSSION 

In this study, we globally analyzed SRSFlO-regulated AS 
events based on deep sequencing of the chicken transcrip- 
tome. We found that SRSF10 is involved in all of the 
common modes of AS, with cassette exons being the 
most frequent targets for SRSF10 regulation. The most 
striking finding of our analysis is that SRSF10 shows 
position-dependent activity. Later, we discuss the 
possible mechanisms by which SRSF10 regulates exon in- 
clusion and exclusion compared with the previously 
reported mechanisms of splicing regulation by other 
RNA regulatory proteins. 

We previously demonstrated that SRSF10 functions as 
a sequence-specific splicing activator by directly binding to 
exonic splicing enhancers. The proposed mechanism for 
SRSF10 activation involves the stabilization of the 
binding of Ul and U2 small nuclear ribonucleoproteins 
(snRNPs) to the 5' splice site and branch site, respectively 
(13). Consistent with this positive regulation, in this study 
we identified that SRSF10 activates a large subset of exon 
inclusions in chicken cell lines, including cassette exons, 
alternative 5' exons and alternative 3' exons. Motif 
analysis and functional studies conducted on a model 
minigene provided strong evidence that when binding 
within an alternative exon, SRSF10 acts as an activator 
to promote exon inclusion. This was further corroborated 
by previous findings in which SRSF10 was shown to 
enhance the inclusion of a triadin alternative 3' exon by 
directly binding to the exon (14). 

However, we also observed that SRSF10 can repress 
exon inclusion in vivo, and thus, sought to determine the 
mechanism underlying this SRSFlO-dependent exon ex- 
clusion. It is possible that when binding occurs within 
flanking constitutive exons, SRSF10 may stabilize the 
interactions of Ul and/or U2 snRNPs with the nearby 
constitutive splice sites, thereby decreasing the competi- 
tiveness of the alternative sites. Kinetic analysis of splice 
site competition has demonstrated that a minor modula- 
tion of splice site recognition might be converted into a 
major functional consequence (30). This suggests that re- 
pression by SRSF10 results from its usual splicing activa- 
tor activity, with partial activation of the flanking 
constitutive sites, preventing the cassette exon from inclu- 
sion. A similar principle has only been previously 
demonstrated for the SR protein SRSF1 (19). It is import- 
ant to note that the manipulation of SRSF10 binding 
sequences could significantly switch an excluded exon 
into an included exon, indicating that SRSFlO-dependent 
activation activity predominates over its repression 
activity. 

Recent genome-wide studies of splicing regulators, 
including the SR proteins, hnRNPs and other RNA regu- 
latory proteins, have revealed that they can all function 
both positively and negatively in regulating exon inclu- 
sion, regardless of their predominant splicing activator 
or repressor activity, which has been previously 



demonstrated using model genes (7-10,31). The opposite 
splicing activity could be attributed to their different 
binding locations in the pre-mRNA, although the pos- 
itional effect exerted by individual splicing regulators 
appears to be fundamentally distinct. However, in some 
cases, the final splicing outcomes are not directly linked to 
the positional effect of individual protein binding in vivo, 
thus emphasizing the combinatory control of AS in vivo. 
Most AS events are likely subjected to regulation by 
multiple different splicing regulators, which may act syn- 
ergistically or antagonistically. In this study, we observed 
the differential effects exerted by SRSF10 in modulating 
exon inclusion. Whether this reflects differences in the 
affinity of SRSF10 for target exons or complex mechan- 
isms involving SRSF10 and other regulatory proteins 
remains to be determined 

Our data also revealed that ~70% of the validated 
splicing variants caused by the loss of SRSF10 have not 
been annotated, which indicates that the current chicken 
transcriptome is incomplete. Importantly, these newly 
identified variants do not result from aberrant splicing 
events, as they are also observed to be normally expressed 
in WT cells, albeit at much lower levels than their corres- 
ponding dominant isoforms. As WT cells are more resist- 
ant to ER stress-induced apoptosis than KO cells, this 
raises the possibility that SRSF10 protects cells from 
stress-induced apoptosis by maintaining the normal 
levels of these major isoforms in cells. By contrast, 
SRSF10 deletion results in isoform switching in these 
events and renders cells susceptible to specific apoptosis 
triggers. This possibility was further demonstrated by the 
observation that KO cells with reconstituted SRSF10 re- 
covered the WT splicing patterns for the majority of 
SRSFlO-verified splicing events and considerably rescued 
the ER stress-related defects. How these novel splice 
variants resulting from the loss of SRSF10 contribute to 
ER stress-induced apoptosis is currently unclear and will 
be an important focus of future studies. 

Together, our findings provide mechanistic insights into 
SRSFlO-regulated exon inclusion and exclusion, which 
depend on its binding location with respect to the 
regulated exons within the pre-mRNA. By regulating 
many stress- and apoptosis-related splicing events, our 
findings demonstrate that SRSF10 is required for cell via- 
bility under conditions of stress. 

SUPPLEMENTARY DATA 

Supplementary Data are available at NAR Online. 
FUNDING 

Ministry of Science and Technology of China 
[2012CB524900 and 2012BAK01B00]; National Natural 
Science Foundation [31170753 and 31070704] and the 
One Hundred Talents Program of the Chinese Academy 
of Sciences. China Postdoctoral Science Foundation 
[2013M541564 to W.W.]. Funding for open access 
charge: National Natural Science Foundation [31170753]. 

Conflict of interest statement. None declared. 



4030 Nucleic Acids Research, 2014, Vol. 42, No. 6 



REFERENCES 

1. Goodell,M.A., Brose,K., Paradis,G., Conner,A.S. and 
Mulligan,R.C. (1996) Isolation and functional properties of 
murine hematopoietic stem cells that are replicating in vivo. 
J. Exp. Med., 183, 1797-1806. 

2. David,C.J., Chen,M., Assanah,M., Canoll,P. and ManleyJ.L. 
(2010) HnRNP proteins controlled by c-Myc deregulate pyruvate 
kinase mRNA splicing in cancer. Nature, 463, 364-368. 

3. Singh,R.K. and Cooper,T.A. (2012) Pre-mRNA splicing in disease 
and therapeutics. Trends Mol. Med., 18, 472-482. 

4. Chen,M. and ManleyJ.L. (2009) Mechanisms of alternative 
splicing regulation: insights from molecular and genomics 
approaches. Nat. Rev. Mol. Cell Biol, 10, 741-754. 

5. Nilsen,T.W. and Graveley,B.R. (2010) Expansion of the 
eukaryotic proteome by alternative splicing. Nature, 463, 457-463. 

6. ChenJ.J., Lin,K.K., Huque,M. and Arani,R.B. (2000) Weighted 
p-value adjustments for animal carcinogenicity trend test. 
Biometrics, 56, 586-592. 

7. Feng,Y., Chen,M. and ManleyJ.L. (2008) Phosphorylation 
switches the general splicing repressor SRp38 to a sequence- 
specific activator. Nat. Struct. Mol. Biol, 15, 1040-1048. 

8. Vethantham,V., Rao,N. and ManleyJ.L. (2008) Sumoylation 
regulates multiple aspects of mammalian poly(A) polymerase 
function. Genes Dev., 22, 499-511. 

9. Llorian,M., Schwartz,S., Clark,T.A., Hollander,D., Tan,L.Y., 
Spellman,R., Gordon,A., Schweitzer,A.C, de la Grange,P., Ast,G. 
et al. (2010) Position-dependent alternative splicing activity 
revealed by global profiling of alternative splicing events regulated 
by PTB. Nat. Struct. Mol. Biol., 17, 1114-1123. 

10. Finck,B.N., Gropler,M.C, Chen,Z., Leone,T.C, Croce,M.A., 
Harris,T.E., LawrenceJ.C. Jr and Kelly,D.P. (2006) Lipin 1 is an 
inducible amplifier of the hepatic PGC-lalpha/PPARalpha 
regulatory pathway. Cell Metab., 4, 199-210. 

11. Biondo,P.D., Brindley,D.N., Sawyer,M.B. and Field,C.J. (2008) 
The potential for treatment with dietary long-chain 
polyunsaturated n-3 fatty acids during chemotherapy. J. Nutr. 
Biochem., 19, 787-796. 

12. Shin,C, Feng,Y. and ManleyJ.L. (2004) Dephosphorylated 
SRp38 acts as a splicing repressor in response to heat shock. 
Nature, 427, 553-558. 

13. Nunez,E., Kwon,Y.S., Hutt,K.R., Hu,Q., Cardamone,M.D., 
Ohgi,K.A., Garcia-Bassets,L, Rose,D.W., Glass,C.K., 
Rosenfeld,M.G. et al. (2008) Nuclear receptor-enhanced 
transcription requires motor- and LSD 1 -dependent gene 
networking in interchromatin granules. Cell, 132, 996-1010. 

14. Feng,Y., Valley,M.T., LazarJ., Yang,A.L., Bronson,R.T., 
Firestein,S., Coetzee,W.A. and ManleyJ.L. (2009) SRp38 
regulates alternative splicing and is required for Ca(2+) handling 
in the embryonic heart. Dev. Cell, 16, 528-538. 

15. Alves,G. and Yu,Y.K. (2011) Combining independent, weighted 
P-values: achieving computational stability by a systematic 
expansion with controllable accuracy. PLoS One, 6, e22647. 

16. Ray,D., Kazan,H., Cook,K.B., Weirauch,M.T., Najafabadi,H.S., 
Li,X., Gueroussov,S., Albu,M., Zheng,H., Yang,A. et al. (2013) 



A compendium of RNA-binding motifs for decoding gene 
regulation. Nature, 499, 172-177. 

17. Bailey,T.L., Boden,M., Buske,F.A., Frith,M., Grant,C.E., 
Clementi,L., RenJ., Li,W.W. and Noble,W.S. (2009) MEME 
SUITE: tools for motif discovery and searching. Nucleic Acids 
Res., 37, W202-W208. 

18. Li,H., Wang,Z., Zhou,X., Cheng,Y., Xie,Z., ManleyJ.L. and 
Feng,Y. (2013) Far upstream element-binding protein 1 and RNA 
secondary structure both mediate second-step splicing repression. 
Proc. Natl Acad. Sci. USA, 110, E2687-E2695. 

19. HanJ., DingJ.H., Byeon,C.W., KimJ.H., Hertel,K.J., Jeong,S. 
and Fu,X.D. (2011) SR proteins induce alternative exon skipping 
through their activities on the flanking constitutive exons. Mol. 
Cell. Biol, 31, 793-802. 

20. Akerfelt,M., Trouillet,D., Mezger,V. and Sistonen,L. (2007) Heat 
shock factors at a crossroad between stress and development. 
Ann. N. Y. Acad. Sci., 1113, 15-27. 

21. Maloy,K.J. and Powrie,F. (2011) Intestinal homeostasis and 
its breakdown in inflammatory bowel disease. Nature, 474, 
298-306. 

22. Porn-Ares,M.L, Samali,A. and Orrenius,S. (1998) Cleavage of the 
calpain inhibitor, calpastatin, during apoptosis. Cell Death Differ., 
5, 1028-1033. 

23. Chan,C.W., Lee,Y.B., UneyJ., Flynn,A., TobiasJ.H. and 
Norman,M. (2007) A novel member of the SAF (scaffold 
attachment factor)-box protein family inhibits gene expression and 
induces apoptosis. Biochem. J., 407, 355-362. 

24. Ladomery,M. (2013) Aberrant Alternative Splicing Is Another 
Hallmark of Cancer. Int. J. Cell Biol, 2013, 463786. 

25. Ferraiuolo,M.A., Basak,S., DostieJ., Murray,E.L., 
Schoenberg,D.R. and Sonenberg,N. (2005) A role for the eIF4E- 
binding protein 4E-T in P-body formation and mRNA decay. 

/. Cell Biol, 170, 913-924. 

26. Penkert,R.R. and Kalejta,R.F. (2012) Tale of a tegument 
transactivator: the past, present and future of human CMV pp71. 
Future Virol, 7, 855-869. 

27. Temme,N. and Tudzynski,P. (2009) Does botrytis cinerea Ignore 
H(2)0(2)-induced oxidative stress during infection? 
Characterization of botrytis activator protein 1 . Mol. Plant 
Microbe Interact., 22, 987-998. 

28. Farhan,H., Weiss,M., Tani,K., Kaufman,R.J. and Hauri,H.P. 
(2008) Adaptation of endoplasmic reticulum exit sites to 
acute and chronic increases in cargo load. EM BO J., 21, 
2043-2054. 

29. Perry,M.E. (2010) The regulation of the p53-mediated stress 
response by MDM2 and MDM4. Cold Spring Harb. Perspect. 
Biol, 2, a000968. 

30. Yu,Y., Maroney,P.A., Denker J.A., Zhang,X.H., Dybkov,0., 
Luhrmann,R., Jankowsky,E., Chasin,L.A. and Nilsen,T.W. (2008) 
Dynamic regulation of alternative splicing by silencers that 
modulate 5' splice site competition. Cell, 135, 1224-1236. 

31. Licatalosi,D.D., Mele,A., FakJ.J., UleJ., Kayikci,M., Chi,S.W., 
Clark,T.A., Schweitzer,A.C, BlumeJ.E., Wang,X. et al. (2008) 
HITS-CLIP yields genome-wide insights into brain alternative 
RNA processing. Nature, 456, 464-469. 



